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{probability  density  function  (PDF!  for  the  given  map  data.  The  accuracy  of  the 
analysis  depends  on  the  available  data,  the  estimators  employed,  and  the 
determination  of  appropriate  PDF's.  The  example  which  will  be  discussed 
consists  of  a  large  number  of  terrain  regions  in  an  area  of  eastern  Massa¬ 
chusetts.  This  area  was  chosen  because  experimental  data  in  the  form  of 
electromagnetic  forward  scattering  from  the  surface  are  available.  Because 
of  the  number  of  cases  involved,  the  complexity  of  the  multivariate  height 
distributions,  and  the  type  of  measurements  available,  only  a  single  observa¬ 
tion  of  the  multivariate  data  is  used  in  the  present  analysis.  The  results  of 
the  statistical  analysis  for  mean  height,  variance,  correlation,  PDF,  and  a 
geologic  feature  characterizing  each  subregion  are  presented.  These  are  used 
in  a  single -roughness -scale  electromagnetic  scattering  form^tedion  to  allow 
comparison  with  the  experimental  data.  A  computer  program  Has  been 
developed  to  determine  the  amount  of  specular  and  diffuse  multipath  power 
reaching  a  monopulse  receiver  from  a  pulsed  beacon  and  the  associated  bore- 
sight  pointing  error.  Terrain  inhomogeneities  and  multiple  specular  reflec¬ 
tion  points  are  included  in  the  analysis.  Results  are  compared  with  data. 
Currently,  the  terrain  analysis  is  being  extended  to  apply  to  the  two-scale 
theory  for  scattering  with  improved  techniques  for  estimation  of  the  statistical 
parameters  of  the  scattering  model. 
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Prediction  of  Electromagnetic  Scattering  for  Rough  Terrain 
Using  Statistical  Parameters  Derived  From 
Digitized  Topographic  Maps 


I.  INTRODUCTION 

The  characteristics  of  electromagnetic  signals  scattered  from  rough  terrain 

1  2 

include  contributions  from  clutter  return  and  multipath  return.  ’  These  two 
aspects  can  be  described  by  the  theory  of  scattering  from  rough  surfaces  if  proper¬ 
ties  of  the  terrain  such  as  probability  density  function  (PDF)  for  the  surface  height 
distribution,  the  covariance  matrix,  R,  the  variance  in  surface  height,  cr  ,  and  the 
complex  dielectric  constant  characterizing  the  surface  are  known.  The  numerous 

theoretical  models  of  EM  wave  scattering  from  rough  surfaces,  for  example, 

1  2  3  4 

Beckmann  and  Spizzichino,  Ruck  et  al,  lx>ng,  '  and  Brown,  all  relate  the 
normalized  cross  section  of  terrain  to  the  foregoing  parameters  characterizing 
the  rough  surface. 

In  this  report,  the  physical  parameters  of  the  rough  surface  are  obtained 
from  digitized  terrain  maps  (furnished  by  the  Electromagnetic  Compatibility 

(Received  for  publication  2  September  1980) 

1.  Beckmann,  P. ,  and  Spizzichino,  A.  ( 1 963)  The  Scattering  of  Electromagnetic 

Waves  from  Rough  Surfaces,  Macmillan  Co. 

2.  Ruck.  Ci.T.,  Barrick,  D.  E.  ,  Stuart,  W.  D.  ,  and  Krichbaum,  C.  K.  (1970) 

Radar  Cross  Section  Handbook,  Vol.  2,  Plenum  Press. 

3.  Bong,  M.W.  (1975)  Radar  Reflectivity  of  Land  and  Sea,  Lexington  Books. 

4.  Brown,  CBS.  (1978)  Baekscattering  from  a  Gaussian  -  distributed  perfectly 

conducting  rough  surface,  IEEE  Trans,  on  Antennas  and  Prop. 

AP-26(No.  3>:472. 
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Analysis  Center,  ECAC,  and  the  Defense  Mapping  Agency,  DMA).  Estimation 
theory  is  employed  to  specify  the  corresponding  statistical  parameters.  A 
hypothesis  testing  procedure  determines  the  probability  density  function  (PDF)  for 
the  surface  heights. 

The  specific  problem  used  as  an  example  is  that  of  characterizing  a  large 

2 

terrain  region  considered  to  be  made  up  of  smaller  subareas  (~  4  km  ).  The  main 
feature  of  interest  is  the  distribution  of  heights  within  these  subregions.  Each 
subregion  is  characterized  by  a  geologic  code  and  several  statistical  parameters. 

In  particular,  we  are  concerned  with  being  able  to  associate  a  PDF  with  the  range 
of  heights  (z.)  in  the  subregions  and  to  determine  parameters  that  make  the  general 
PDF  explicit.  The  data  elements  z.  =  z.lx^,  Y f,  )»  where  i  =  1,  2,  3, . . .  N,  and  N 
is  the  total  number  of  grid  points  in  the  x  -  y  plane  constituting  the  subregion. 

Here,  x.  denotes  the  kth  equally  spaced  x-value  along  the  x-axis  and  y .  denotes 

th  ^  r"*  £  i*“ 

the  £  v -value  along  the  y -axis,  where  k  =  1,2,..,  \N  and  ] C  =  1,2,,,.  \N, 

Thus,  the  N  points  are  distributed  in  the  x-y  plane  in  order  to  form  a  rectangular 
grid.  The  covariance  matrix  can  be  assumed  to  have  the  form: 

R  a2  exp  (-  t"  / T2) 

where  T  is  the  correlation  length  and  for  this  class  of  data  sets,  the  form  of  t 

mn 

is 


T  =  (x  -x  )  +  (y 

mn  m  n  m 


V 


The  motivation  for  assuming  a  covariance  matrix  of  this  form  is  that  it  leads  to  a 
tractable  mathematical  expression  for  the  incoherent  power  scattered  when  an 
electromagnetic  wave  is  reflected  from  a  rough  surface.  Similarly,  in  the  hypoth¬ 
esis  testing  procedure,  the  binary  decision  process  involves  the  case  where  the 
PDF  is  either  an  N -variate  Gaussian  or  exponential.  This  specialization  is  also 

1  2 

motivated  by  the  theory  of  electromagnetic  wave  scattering  from  rough  surfaces.  ’ 

The  normalized  radar  cross  section  of  the  rough  surface  used  in  this  report 
6  7  8 

was  derived  by  flagfors,  Barrick,  and  Semenov.  This  cross  section,  Oq,  is 


5.  Whalen,  A.  D.  (197  1)  Detection  of  Signals  in  Noise,  Academic  Press. 

6.  Hagfors,  T.  ( 1964)  Backscattering  from  an  undulating  surface  with  applications 

to  radar  returns  from  the  moon,  J.  Geophys.  Res.  69:377  9. 

7.  Barrick,  D.  E.  (1968)  Relationship  between  slope  probability  density  function 

and  the  physical  optics  integral  in  rough  surface  scattering,  Proc.  IEEE 
56:1728.  - 

8.  Semenov,  B.  ( 1965)  Scattering  of  electromagnetic  waves  from  restricted 

portions  of  rough  surfaces  with  finite  conductivity,  Radiotekh.  i  Elektron, 
10:1952.  - 

/VK\ 
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incorporated  into  a  computer  program  which  calculates  the  amount  of  specular 
and  diffuse  multipath  power  entering  a  monopulse  receiver  from  a  beacon  located 
over  rough  terrain  (see  Figure  1).  The  computer  program  also  calculates  ihe 
error  in  boresight  pointing  accuracy  of  a  monopulse  receiving  antenna  due  to 
noise  and  diffuse  multipath.  Among  other  things,  the  computer  program  takes  into 
consideration,  (1)  the  spatial  nonuniformity  of  the  rough  earth  (that  is,  the  pre¬ 
ceding  characterization  parameters),  (2)  nonuniformities  in  the  glistening  surface, 
(3)  finite  pulse  length  of  the  beacon,  (4)  antenna  elevation  power  pattern  of  the 
monopulse  receiver,  (5)  multiple  specular  reflection  points  due  to  unevenness  in 
surface  height,  (6)  interference  between  direct  signal  and  multiple  specularly 
reflected  signals,  and  (7)  finite  azimuthal  beamwidth  of  transmitter  and  receiver. 

9 

Finally,  for  the  case  of  normally  distributed  surface  heights,  Sancer's  results 
are  employed  to  describe  the  effects  of  shadowing  and  Brown's10  general  back- 
scatter  shadowing  results  have  been  extended  to  handle  the  exponential  case. 


PULSEO  RADAR  BEACON 


DIFFUSE 
MULTIPATH  - 


DIRECT  SIGNAL - 


SPECULAR  MULTIPATH  • 


MONOPULSE  RECEIVER 


ROUGH  TERRAIN 

Figure  1.  Reflection  of  Radar  Waves  From  Rough  Terrain 

The  results  of  the  program  describe  the  effect  of  the  terrain  on  the  electro¬ 
magnetic  signal.  The  final  outputs  include  total  coherent  and  diffuse  power  levels, 
and  the  induced  boresight  error.  The  data  output  from  the  computer  program  for 
the  sum  pattern  coherent  power  and  angular  error  in  boresight  is  compared  with 
experimental  data  taken  by  personnel  at  the  MIT  Lincoln  Laboratory  Discrete 
Address  Beacon  System  (DABS)  test  site.  11 

9.  Sancer,  M.  I.  (1969)  Shadow -corrected  electromagnetic  scattering  from  a 

randomly  rough  surface,  IEEE  Trans,  on  Antennas  and  Prop.  _AJP^17  ;577 -585 

10.  Brown,  G.S.  (1980)  Shadowing  by  non-Gaussian  random  surfaces.  Proceedings 
of  the  Second  Workshop  on  Terrain  and  Sea  Scatter,  George  Washington 


university,  Washington,  u.  u. 

11.  McGarty,  T.  P.  (1975)  The  Statistical  Characteristics  of  Diffuse  Multipath 


and  its  Effect  on  Antenna  Performance,  AD- 


2.  STATISTICAL  ANALYSIS 


In  our  statistical  analysis  of  terrain  heights,  we  carry  out  a  series  of  opera¬ 
tions:  First,  we  propose  some  multivariate  probability  density  functions  that  may 
represent  the  observed  distribution  of  height  values;  next  we  use  the  available  data  to 
generate  appropriate  estimators  of  the  parameters  of  the  respective  densities; 
finally,  we  conduct  a  hypothesis  test  to  ascertain  which  density  function  is  more 
likely  to  have  produced  the  observed  height  data. 

The  two  PDF's  for  the  heights  are  multivariate  Gaussian  or  exponential.  The 

12 

Gaussian  has  a  well-known  form: 


PG<Z1’  z2* 


SN)  =  [ 


(2t,) 


N/2 


[1  »  T  - 1 

-  (z  -  p)  R 


where  R  represents  the  covariance  matrix.  For  our  case,  we  assume  equal  means 
~  2  2 
(p),  variances  (a  )  and  that  the  correlation  function  (pj.  =  a  /a  )  has  a  Gaussian 

dependence  on  the  separation  between  points.  The  next  aspect  is  to  develop  a 

similar  form  for  the  exponential: 


PE<ZrZ2’ 


. . . , )  =  C  j  exp 


<z' 


The  two  coefficients  (Cj,  C2>  have  to  be  determined.  This  will  result  in  a  form 
that  satisfies  the  requirements  for  a  PDF.  To  evaluate  the  coefficients,  we  use 
the  properties  that  the  zeroth  moment  integral  of  a  PDF  is  equal  to  unity  and  the 
second  moment  integral  is  equal  to  the  variance.  We  thus  obtain  a  form  for  the 
multivariate  exponential: 


PE  (zl'  Z2' 


V 


[•v 


(277) 


N  -  1 


(N+  1)N,/2  exp 


ft 


<N  +  1)  (z  -  p>T  R"1  (z 


«>]  '"j 


In  order  to  decide  which  of  these  two  PDF's  is  more  appropriate  for  the  data, 
we  must  next  establish  estimators  for  the  parameters  of  the  densities:  the  means, 
variance,  and  covariances.  The  complexity  of  correlated  multivariate  analysis  and 
the  computational  limitations  associated  with  the  available  data  formats  caused  us 
to  select  estimators  that  have  an  intuitive  appeal  based  on  their  form  rather  than  a 
rigorous  derivation.  From  our  assumption  of  equal  means  and  variances,  we  use 
the  sample  mean  as  the  estimator  for  the  mean  height  and  the  sample  variance  as 
the  variance  estimate: 


12.  Mood,  A.  M.  ,  and  Graybill,  F.  A.  (1963)  Introduction  to  the  Theory  of  Statistics, 
McGraw-Hill. 


«  • 


VT  VFT 

r=(\in)  5Z  z(xi’  ■yjl  (mean)  and 

i=  1  n=  1 


a  =  (l/N) 


yjN  Vn 

jj)  ^  ^zi  ‘  (variance>- 


i=  1  2=  1 


The  procedure  for  the  covariance  matrix  estimator  is  more  complicated.  A 
correlation  length,  T,  is  defined  as  the  separation  at  which  a  normalized  covariance 


function  Cz?  has  decreased  to  the  value  e 


where  C 


■  y  /<x  ,  and  y  is  the 
'  zz  '  zz 


estimator  for  the  covariance.  The  data  are  used  to  determine  the  estimate  of  T  in 
this  fashion  and  then  the  complete  covariance  estimator  is  formed  from  the  relation: 

>zztni'n)  =  52exp(-Tmn/T2) 

(separate  X  and  V  relations  are  calculated).  The  form  of  the  covariance  (X  direc¬ 
tion)  used  to  determine  T  is 


C  00  =  (1/Nff  ) 
zz 


VN  VN-k 

E  E 

- 1  1=1 


Z ,  .  Zi .  i  —  z 

i.l  i  +  kj 


VN  VN-k 

E  E 

.1=1  i=J 


Z  - 

Vn- 


where  k  represents  the  separation  distance.  A  least  squares  fit  of  the  C'zz  values 
to  a  parabola  is  then  used  to  find  T.  Under  the  assumptions  that  have  been  made, 
the  above  estimators  are  similar  to  those  found  in  Jenkins  and  Watts.  After 
obtaining  the  values  for  T  in  the  X  and  V  directions,  the  two  values  are  averaged 
to  find  the  final  estimate  of  T  in  a  given  subregion. 

We  now  have  the  two  PDF  forms  and  the  parameters  required  for  them.  The 
next  aspect  is  to  decide  which  PDF  is  the  more  appropriate  for  the  given  data.  One 
final  comment  on  the  parameters  should  be  made.  In  order  to  satisfy  the  restric¬ 
tions  of  the  PDF’s,  it  is  necessary  to  show  that  the  quadratic  form  appearing  in 
both  cases  is  positive  definite.  This  has  been  demonstrated  for  the  above  cases  by 

making  use  of  the  Gaussian  form  assigned  to  the  covariance  matrix  elements. 

14 

Details  can  be  found  in  the  report  by  Lennon  and  Papa. 

The  form  of  the  hypothesis  test  used  here  is  based  on  the  maximum  a  posteriori 
probability  criterion.  This  is  equivalent  to  a  minimum  error  probability  criterion. 
We  assign  hypothesis  Hj  to  the  Gaussian  case  and  hypothesis  Hq  to  the  exponential. 
Then  the  likelihood  ratio  parameter, 

13.  Jenkins,  G.  .VI. ,  and  Watts,  D.G.  (1968)  Spectral  Analysis  and  its  Applications, 

Holden -Day. 

14.  Lennon,  J.  F. ,  and  Papa,  R.  J.  (1980)  Statistical  Characterization  of  Rough 

Terrain,  RADC-TR-80-9,  RADC/EE  Hanscom  AFfe,  Massachusetts, 
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Let  P  (Hq)  be  the  probability  that  hypothesis  Hq  is  true.  Then  the  decision  rule 
may  be  written  as:  Choose  Hj  if 


P(H  ) 
o 


For  our  case,  we  assume  that  it  is  equally  probable  that  hypothesis  Hj  orH^  is 
true  and  the  decision  rule  reduces  to  whether  or  not  A  2  1.  Note  that  it  may  be 
possible  to  alter  the  probability  that  H  is  true  based  on  external  evidence  (such 
as  the  type  of  terrain). 

When  the  specific  forms  for  the  two  PDF's  are  introduced  into  this  relation  it 


=  (z  -  p)1  H  1  (z  -  p)j 


For  convenience  we  rewrite  the  test  in  logarithmic  form  and  obtain  the  result  that 
Hj  is  true  if 


or 

Q  2  VN+  1  -  j^2  ln[r(a^-!)j  T  (\'  +  D-  In;  -  \ 

For  the  actual  cases  N  =  100  and  the  specific  result  is  that  we  decide  that  the 
terrain  heights  in  a  given  region  are  from  a  Gaussian  PDF  if 

85.01  '  [<z  -  /I)  r  if  1  </.  -  ji  )j  S  U8.37 

■>  ■> 

and  conversely  the  points  are  exponential  if  Q“  >  1  18.  J7  or  <  85.  01. 

These  procedures  have  been  applied  to  a  specific  site  to  allow  comparison  of 
theoretical  electromagnetic  calculations  based  on  the  terrain  eha factorization  with 
actual  experimental  data  available  for  the  territory  described. 
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1  If.  Te  :  •:  r« » :•  tin  cha  r.ieter;/  die;  in  <  ,.-tert:  V.  is.,  ,'i  < -ft s. 

I.’  ■  1. 1  s  •  •.  .ip!  .  ■•iiit  ■  it-;  on  eh  it  rom  .  ._=  sift  i  ■  s  .  t  *  •  -  f  i  r*i?  ?■•<>:  ■  i  t  :i;<1  terrain, 

V  ••  t  •  ■  -o  .  oil. pa  re  t.  ■ -■<  i*. u:  il  imI.-u!  ’  i  ns  •><  it  :i  .-xpi  "  i  n :  •  •  nl .  i  t  ro- 

-  \  i!  r  r-.-j  round  this  Di.e..-  rote  Mr  li-  ..-un  svst.  :::  IJABM  site 

->  lii.-.-iitn  i  is  :•:•••  ,  w-s  4  km  long  and  42.  d  ki;  wide.  It  u,  s  subdivided 

-•  •  si;.  >  1 1  ■  ■  ri  rt.  lit:. ilar  ueils,  each  with  sides  of  20"0  ::  bv  1  o  a  2  \v  ,  t  ; ,  ,  Mcthe!' 

S-.  l'r  ii  .  ;  sii  >:i  if •  •  ;  e.h  cel  1  into  a  i  0  -  1  0  g  rid  of  points.  \  : 1 1  a  ha  -  c  .  ‘up-  ig  r  uphu 

1  ;i a  ,  a.t  ;  ,>.l-  *  \  |  ,.  *  >1  j  _  .  |  ea  is  available  Jt  till'  i  1 1  ■  i'T  !'■  ■  1 1 ..  -  r  1 1  •  •*  i  (  (  I  'HI  p  a  I  i  hi !  it  v  V;.  I !  VS  is 
crd-  :  i  l.'C  »('>,  prep.: I'i'u'  from  Defense  M  .pping  \getics  dill  \t  supplied  digitized 
Iff  ruin  maps  at  1:  2S0, 000  scale  size. 

riv  statistical  lata  for  e. .eh  coll  has  been  re.  .run!  .  -n  .  computer  tape  fur  iso 
with  the  program  for  the  electromagnetic  analysis,  ilach  cell  is  rcprt  seated  bv 
seven  descriptors.  The  first  two  entries  are  the  ■  X ,  \  I  coordinates  for  the  renter 
of  the  eel]  I  tin  •  origin  of  tin-  coordinate  system  is  taken  as  the  center  of  tin-  extreme 
southwestern  corner  of  to.  rectangular  region).  The  next  item  is  the  geological 
cod'  tor  the  c  II.  Tie  predominant  feature  is  woods;  there  are  a  number  of  cells 
containing  elr.-n  "s  of  likes  and  ponds  and  a  few  town  sites  with  associated  cleared 
are  ,s.  This  is  followed  bv  the  mean  nd  variance  of  the  heights  in  the  cell  and  the 
estimated  correlation  l<  ngth  I’.  idle  units  of  length  are  in  meters.  Tin-  final 
quantity  the  result  •>!  the  hypothesis  test.  This  result,  is  presented  in  a  format 
s-.t.-'h  that  the  ti  ight  -•  •  region  are  (l.ms.si.in  when  26  *  TKST  «.  0.  W  hen 

I  I  Si  ;  0  and  the  v.r.  eiee  i  :  very  small,  the  region  is  essentially  a  smooth  surface 
'  ■■(•!  roughness). 

\\  in  n  the  over-  >11  r<  -alts  for  the  region  are  examined.  one  observes  that 
v.!'.  -n  the  magnitude  of  Tl.ST  is  very  large,  correlation  length  is  also  large,  and  for 
those  cases,  1'  i>  comparable  to  one -ha  If  the  cell  size  or  even  larger.  When  this 
occurs,  tin'  determinant  of  tnc  eovariaricc  matrix  K  becomes  very  small.  As  a  re¬ 
sult,  it  becomes,  more  diflu  ult  to  obtain  an  accurate  inverse  of  It  due  to  rapid  build- 

up  of  round-oil  -  rror.  A  rce-ted  difficulty  is  that  for  some  cases  TliST  <  -  1  111.  27. 

) 

Tins  would  be  pussibb  uni .  if  the  quadratic  form  Ct-~  is  not  positive  definite.  Thus, 
further  machine-induced  errors  arc  present,  and  those  results  can  not  be  considered 
valid.  I  urthor  consideration  is  being  given  to  the  statistical  aspects. 

in  order  to  use  these  results  in  the  rough  surface  electromagnetic  calculations, 
one  additional  aspect  should  be  noted,  l  or  the  types  of  geological  features  that 

describe  the  respective  regions,  data  exists  on  the  associated  complex  dielectric 

If),  2 

constants  at  microwave  1  requeneies. 
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4.  ELECTROMAGNETIC  SCATTERING  ANALYSIS 


The  radar  cross  section  of  terrain  is  normalized  with  respect  to  the  average 

area  illuminated  by  the  radar.  A,  The  normalized  cross  section  O  may  be  divided 

1 6  ^ 

into  three  general  categories:  (1)  The  slightly  rough  surface;  (2)  the  very  rough 

1678  17  18  4 

surface;  '  '  ’  and  (3)  the  multiple  scale  rough  surface.  ’  *  This  report  is 

concerned  with  the  second  category,  that  is,  the  very  rough  surface. 

2 

Ruck  et  al  give  expressions  for  the  average  bistatic  rough  surface  cross 

section  a  under  the  following  four  assumptions:  (1)  the  radius  of  curvature  of  the 

surface  irregularities  is  larger  than  a  wavelength;  (2)  the  roughness  is  isotropic  in 

both  surface  dimensions;  (3)  the  correlation  length  is  smaller  than  either  the  X 

or  Y  dimension  of  the  sample  subregion;  and  (4)  multiple  scattering  is  neglected. 

2 

Using  the  notation  of  Ruck  et  al  one  notes  that  the  expression  for  cr  becomes: 


a  = 
o 


pq 


2 


J 


whe  re 


J  = 


„  exp 


a  « 


for  a  Gaussian  bivariate  surface  height  probability  density  function  and 


for  an  exponential  surface  height  probability  density  function.  The  scattering 
matrix  elements  8  are  given  by 

pq 


( 1  +  cos  2of )  K ||  (a) 

8  =  (cos'g,  T  cos  FT  (vertical  polarization) 

l  s 


16.  Peake,  W.  H.  (1959)  The  interaction  of  Electromagnetic  Waves  with  Some 

Natural  Surfaces,  Antenna  Laboratory,  Ohio  State  University  Report 
- 

17.  Wright,  J.W.  (1968)  A  new  model  for  sea  clutter,  IEEE  Trans,  on  Antennas. 

and  Prop.  AP-16: 2 1 7-223. 

■  ■  ■  ■ 

18.  Kuks,  1.  (1966)  Contributation  to  the  theory  of  radio  wave  scattering  on  the 

perturbed  sea  surface,  lz.  Vyssh,  Ucheb.  Zaved.  Radiofiz  5:876. 
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(horizontal  polarization) 


1 


^  hh  = 

R  ||  (a) 

R_L  (a) 

where 

9.  -  angle  of  incidence  (with  respect  to  surface  normal) 
0^  =  angle  of  scattering  (with  respect  to  surface  normal) 
!x  =  sin  0.  -  sin  0g,  =0.0.  Iz  =  ~  c°s  -  1.0. 


(1  +  cos  2a)  to) 


(cos  9.  +  cos  9  ) 

1  s 

e  cos  a  - 

fe  -  sin*  a 
r 

€  cos  a  +  »J 

i - r-  2  — 

r €  -  sin  a 

r  Y 

r 

cos  a  — 

- “2 — 

-  sin  a 

- ;  -  ^-2- 

cos  a  +•  */e  -  sin  a 


r 


and 


a 


-( 


e.  +e 

1  s 
2 — 


)• 


Here,  e  is  the  relative  complex  dielectric  constant  of  the  surface,  the  subscript 
r 

||  refers  to  the  E -field  in  the  plane  of  incidence,  and  the  subscript  j_  refers  to 
the  E -field  normal  to  the  plane  of  incidence.  These  simplified  forms  of  Ruck's 
expression  follow  from  the  assumption  that  the  receiver  is  far  from  the  trans¬ 
mitter  so  that  the  portion  of  the  "glistening  surface"  that  contributes  to  the  diffuse 
multipath  is  a  long,  narrow  strip  extending  between  the  transmitter  and  receiver. 
This  assumption  allows  us  to  make  the  approximation  that  the  azimuthal  scattering 
angle  <pg  ~  0.  0.  To  describe  the  effect  of  shadowing,  the  previous  expression  for 
a  is  multiplied  by  the  appropriate  shadowing  function. 

It  was  explained  in  the  Introduction  how  a  computer  program  was  developed  to 

incorporate  the  expression  for  a  into  an  integral  over  the  glistening  surface.  In 

o 

order  to  calculate  the  amount  of  coherent  power  (specular  plus  direct)  and  diffuse 
multipath  power  reaching  a  monopulse  receiver  from  a  beacon  located  over  rough 
terrain  (see  Figure  1),  the  computer  program  uses  the  previously  described  data 
tape  of  the  statistical  parameters  for  a  particular  site  as  an  input.  Other  input 
variables  characterize  the  transmitter,  the  environmental  aspects,  and  the  re¬ 
ceiver.  Transmitter  values  include  the  gain  of  the  transmitter,  the  polarization  of 
the  transmitted  wave,  the  peak  power  of  the  transmitted  pulse,  the  pulse  length, 
and  the  wavelength  of  the  signal.  External  inputs  include  the  complex  dielectric 


constant  of  each  type  of  geological  region,  the  coordinates  (I,."- 1 : u» l - ■  arc  Ion. -male  • 
of  the  rr.onopulse  receiver,  the  initial  an.;  anal  |>osr. ton  •>:  the  ur. •••.,•;  y 

the  transmitter,  a  parameter  to  control  the  el'i'o  of  su.viow i n /,  .a,  :  .  ■  n  - 

of  the  aircraft.  Receiver  data  are  the  height  oi  the  ren-iv  r,  Ice  : .  -lie 

receiver,  the  front-end  receiver -noise  figur-,  the  antenna  gait.  ; -  .-  cr  inn 

difference  patterns  of  the  monopulse  receiver,  its  azimuthal  li'-anv.  mi r.,  •  ne 
sampling  frequency  of  the  receiver,  the  transmission  line  loss  factor  of  the  it/'c-. 
connecting  ttie  antenna  to  receiver,  and  the  difference  pattern  slope  near  she  born-  - 
sight  axis. 

From  a  knowledge  of  the  initial  and  final  positions  of  the  aircraft  and  the  air¬ 
craft  speed,  the  computer  program  first  calculates  the  trajectory,  as  a  function 
of  time.  Then,  from  a  modified  form  of  the  radar  range  equation,  the  electric 
field  intensity  of  the  direct  signal  at  the  receiver  is  calculated  at  fixed  time  inter¬ 
vals  (sampling  time  is  an  input  variable  to  the  program).  For  the  specular  multi¬ 
path  rays,  all  possible  specular  reflection  points  between  the  transmitter  and 
receiver  are  determined  for  each  position  of  the  transmitter.  Multiple  specular 
reflection  points  due  to  unevenness  in  surface  height  and  possible  rav  blocking  are 
taken  into  consideration.  The  finite  dielectric  constant  of  the  eartli  at  each  specular 
point,  the  antenna  eievatjon  pattern  (receive),  and  the  surface  roughness  are  also 
accounted  for  in  calculating  the  phase  and  amplitude  of  each  specular  multipath  rav. 
At  each  point  on  the  transmitter's  trajectory,  the  total  coherent  power  for  the  sum 
and  difference  channels  in  the  monopuiso  receiver  is  calculated.  Various  assump¬ 
tions  on  the  roughness  effects  have  been  examined.  Those  are  discussed  in 
Section  5. 

The  diffuse  multipath  power,  Pj  entering  the  receiver  is  obtained  from  tin- 
equation 
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In  Hit-  tiuinerii  al  interrulion  of  (In-  equation  lor  i'|(.  a.  i-ount  is  l.iiim, 
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azimuthal  beamwidth. 

1’  ■  signal  to  interference  ratio  in  the 

difference  i  hairnet  . 

i  iihfmu  (tower  \n  sum  channel, 

normalized  monopulse  slope  (obtainable  from  sum 
and  difference  patterns). 


The  output  of  the  computer  program  consists  of  azimuthal  angular  error  in 
boresight  due  to  noise,  the  error  due  to  noise  plus  diffuse  multipath,  the  total  sum 
pattern  coherent  power,  the  total  diffuse  power,  signal -to-noise  ratio,  signal -to- 
interference  ratio,  and  range  from  transmitter  to  receiver.  In  the  conclusions, 
the  computer  output  for  the  analysis  of  a  particular  site  is  compared  with  some 
experimental  data. 


5.  RI.SI  |.TS  VNI)  COM.U  SIONS 

Two  theoretical  results,  the  coherent  power  in  the  sum  channel  of  a  mono¬ 
pulse  antenna,  1  *  a,lf*  •he  standard  deviation  in  boresight  pointing  aecuraev, 

<jg,  arc  compared  with  the  experimental  data  of  McGarty,  **  taken  at  the  Discrete 
Address  Beacon  System  (DABS)  site. 

The  DABS  receiver  was  an  L-band  rotatable  array.  The  beacon  was  located 
on  a  l  -10  aircraft  which  flew  a  number  of  radial  trajectories  mward  and  away 
from  the  monopulse  receiver.  Data  were  recorded  for  abou1  100  flights  for  differ¬ 
ent  aircraft  heights,  different  radial  flight  trajectories,  and  different  receiver- 
antenna  tilt  angles.  The  conditions  under  which  the  data  were  taken  are  listed  in 
Table  1. 

In  order  to  make  the  comparisons,  various  manipulations  were  required.  In 
the  equation  for  a ,  the  normalized  slope  in  the  difference  pattern  near  boresight 
k  was  obtained  from  graphs  of  the  sum  and  difference  patterns  of  the  monopulse 
antenna.  Within  the  accuracy  of  these  graphs,  it  appears  that  1.  5  -  -  1.7. 

Also,  instead  of  integrating  the  difference  pattern  in  azimuth  over  the  glistening 
surface  to  obtain  the  diffuse  multipath  power  received,  an  average  value  ( 18.  5  dB) 
for  the  difference  channel  gain  was  used  over  a  3*  beamwidth. 


6 


Table  1.  Experimental  Conditions  for  Discrete  Address 
Beacon  Svstem  Tests 


Front  end  receiver  noise  figure 

2  >IB 

(lain  of  monopulse  receiver  (sum  pattern) 

22.  5  d B 

(lain  of  transmitter 

4  dB 

Height  of  receiver 

101  m 

Height  of  transmitter 

1220  m 

Signal  polarization 

vertical 

Beak  transmitter  power 

250  W 

Pulse  length 

20  psee 

Azimuthal  beamwidth  (receiver) 

2  5 

Wavelength 

0.  275  m 

Transmission  line  loss  factor 

2  dB 

Tlie  experimental  results  are  illustrated  in  Figures  2  and  2.  Figure  2  is  a 
plot  of  the  sum  signal  (l^  ^,)  versus  range  of  transmitter  to  receiver  in  nmi.  Fi 
ure  2  is  a  plot  of  azimuth  error  (ct^>  versus  range  in  nmi. 


lb  70 
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Figure  2.  Experimental  Data;  Sum  Signal  vs  Range 
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RANGf.  intro  i 

Figure  3.  Experimental  Data:  Azimuth  Error  (ct  )  vs  Hange 


When  the  variance  associated  with  the  distribution  oi  heights  over  the  whole 
region  is  used,  the  coherent  power  from  the  specular  ravs  is  minimal  and  die  i  ,• 
suit  is  just  the  expected  direct  rav  fall -off  of  ixiwer  with  range.  I  his  is  consistent 
with  the  scale  of  the  actual  contributing  region  ion  tnc  order  >1  the  first  Fresnel 
zone)  and  the  corresponding  assumptions  it  -mall  .  ir relation  length  in  'lie  s.  alter 
ing  formulation.  The  use  of  more  localized  variani  es  -is so.  i.ate  I  with  the  I  restiei 
zone  dimensions  produces  more  multipath  effects.  Figure  1  shows  the  almlste.i 
power  for  the  sum -signal,  together  with  the  actual  data.  I  he  specular  .mu  intui¬ 
tions  come  from  three  areas  and  the  effect  ot  the  different  to.  al  variances  .i.-signc.i 
to  each  can  be  seen  in  the  changing  behavior  ot  the  signal  along  die  trujectoi  ,• .  |i 
should  he  noted  that  the  theoretical  results  tend  in  lie  iliout  ti  ill  higher  on  the 
average,  probably  due  to  the  fact  that  the  processing  losses  uv  iimkiiowii  an  I  heme 
have  been  neglected.  Figure  5  represents  an  utteinpi  to  correct  mi  dus  i  |.v 

introducing  a  total  processing  loss  of  -a  dll. 

The  next  figures  show  the  results  for  calculated  bores iglit  error  when  both 
extremes  for  k  are  used  and  different  assumptions  made  as  in  shadow  ing  an  I  t  lie 
PDF  of  the  terrain  heights.  For  the  particular  trajectory  of  the  experiment,  ihe 
statistical  analysis  indicates  that  in  all  the  terrain  regions  contributing  to  the 
calculations,  the  heights  are  best  described  as  being  from  exponential  listi  ilmtions. 
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Figure  4.  Sum  Signal  (P  .  ) 
vs  Range:  Theoretical  co 
Results  Without  Processing 

Losses  ( - )  and  Experimental 

Data  ( — ) 


Figure  5.  Sum  Signal  (P  .  ) 
vs  Range:  Theoretical  co 
Results  With  Typical 
Processing  Losses  (---) 
and  Experimental 
Data  ( — ) 


In  Figure  6,  km=  1,  5  and  the  results  for  shadowing  are  shown  for  the  case 

where  processing  losses  are  included.  Comparison  of  Figure  6  with  Figure  3 

shows  that  the  theoretical  boresight  errors  a  &  are  less  than  the  experimental  values. 

This  mav  be  attributable  to  not  integrating  the  true  monopulse  difference  pattern  in 

azimuth  over  the  entire  glistening  surface.  Our  use  of  an  average  18.  5  dB  gain 

over  a  3*  beamwidth  may  have  resulted  in  less  received  diffuse  multipath  power. 

Two  factors  neglected  in  the  program  can  lead  to  increased  boresight  error.  These 

are  the  azimuthal  variation  in  and  the  fine  scale  surface  roughness  contribution. 

However,  even  the  present  agreement  between  theory  and  experiment  is  reasonably 

good.  It  may  be  noted  that  shadowing  reduces  the  error  a ^  only  a  very  small 

amount  for  the  system  parameters  of  this  report. 

Next  we  show  the  results  for  -  1.7.  Figure  7  shows  the  boresight  error  for 

an  exponential  surface  including  the  effect  of  shadowing.  Comparison  with  Figure  3 

shows  that  the  higher  value  for  results  in  slightly  poorer  agreement  between 

theory  and  experiment.  The  conclusion  is  that  for  the  range  of  k  values  used 

m 

here,  there  is  no  significant  effect  on  d  .  To  complete  the  analysis  of  the  effect 
of  PDF,  we  have  the  case  where  all  the  subregions  are  assumed  to  have  Gaussian 
distributed  surface  heights  with  -  1.5.  These  results  are  shown  in  Figure  8. 
Comparison  with  Figure  6  indicates  that  the  assumption  that  the  surface  heights 
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are  all  normally  distributed,  results  in  less  diffuse  multipath  power  entering  the 
receiver.  This  represents  poorer  agreement  with  the  data. 


Figure  6.  Theoretical 
Calculations:  CT g  vs  Range: 
k  =  1,5,  Shadowing 
(Actual  Exponential  Heights) 


Figure  7.  Theoretical 
Calculations:  Oq  vs  Range: 
km  =  1.7,  Shadowing 
(Actual  Exponential  Heights) 


Figure  8.  Theoretical  Calcu¬ 
lations:  Qq  vs  Range:  km  =  1.5, 
Shadowing  (Assumed  Gaussian 
Heights) 


McGartv  also  made  some  comparisons  between  the  experimental  data  and 

theoretical  calculations  based  on  the  rough  surface  scattering  model  proposed  by 
20 

Barriek.  McCarty's  work  involved  the  use  of  the  channel  spread  function,  which 


20.  Barriek,  1).  E.  (1968)  Rough  surface  scattering  based  on  the  specular  point 
theory,  IEEE  Trans,  on  Antennas  and  Frop.  AH-  16:449-454. 
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Lakes  into  consideration  the  fact  that  the  diffuse  multipath  power  has  a  wavenumber 
spectrum  associated  with  it.  His  results  show  reasonably  good  agreement  between 
theory  and  experiment,  but  there  is  no  discussion  of  how  the  surface  correlation 
length  and  the  variance  in  surface  height  were  estimated. 

The  results  presented  in  this  report  are  being  extended.  The  average  power 
restriction  is  being  relaxed.  Better  techniques  for  statistical  parameter  estima¬ 
tion  of  terrain  features  and  alternative  methods  for  the  hypothesis  testing  proce¬ 
dure  are  being  pursued.  Finally,  the  effects  of  mean  surface  tilt  for  each  sub- 
region  and  theoretical  models  of  a  .  which  contain  multiple  scales  of  surface  rough- 

C7 

ness  will  be  included  in  the  scattering  model. 
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